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Abstract 

The work is concerned with the trade-offs between the dimension and the time 
and space complexity of computations on nondeterministic cellular automata. We 
assume that the space complexity is the diameter of area in space involved in 
computation. 

It is proved, that 

1) . Every NCA A of dimension r, computing a predicate P with time complexity 
T(n) and space complexity S(n) can be simulated by r-dimensional NCA with time 
and space complexity 0(T 7 + T S'' : T r ) and by r + 1-dimensional NCA with time and 
space complexity O^T 1 ^ 2 + S), where T and S are functions, constructible in time. 

2) For any predicate P and integer r > 1 if A is a fastest r-dimensional NCA 
computing P with time complexity T(n) and space complexity S(n) , then T = 
O(S). 

3) . If T r _p is time complexity of a fastest r-dimensional NCA computing predi- 
cate P then 

T r+1 , P = 0((T I .,p) 1 -'-/( r + 1 ) 2 ), 
T r _ liP = 0((T r , P ) 1 + 2 /'-). 
Similar problems for deterministic CA are discussed. 
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1. Introduction 

It is well known that nondeterministic computations are more powerful than 
deterministic ones. The interrelation between deterministic and nondeterministic 
time complexity of computations was established by S.Cook in the work [2] where 
he showed the existance of NP-complete problems. However, it is still unknown can 
nondeterministic computations be fulfilled on physical devices or not. In this paper 
we show how the complexity of computations on a nondeterministic device depend 
on its dimension. Note that the similar problem for deterministic computers is open 
(look at the section 6). 

Cellular automata (CA) provide a convenient framework for studies on this prob- 
lem. A cellular automaton is a dynamical system with local interactions operating 
in discrete space and time, and simultaneously CA may be used as a general model 
of computational device. CA were introduced by S.Ulam in the work [7] and J. von 
Neumann in the work [5] and since then various problems pertaining to CA were 
treated in a great many works (look, for example, at [1], [3], [6], [9], [11], [12]). 

Let r > 1 be an integer, Z r be the space, ut be a finite alphabet for possible 
states of any cell i e U . Cellular automaton of dimension r in alphabet w is a 
function of the form: A : u 2r+1 — ► u> . 

A cellular automaton determines the special class of evolutions in Z r so that the 
states of all cells i £ Z r evolve synchronously in discrete time steps according to 
the states of their nearest neighbours. 

It must be mentioned, that a more general approach can be of interest for appli- 
cations, where the function A depends on i or t. It is not our intention to regard 
such possibilities here. 

A configuration is an ensemble of states of all cells at some instant of time. It is 
apparent that an evolution of CA is uniquely determined by the initial configuration. 

If we consider a multifunction instead of A , we obtain the definition of a nonde- 
terministic cellular automaton (NCA). Generally speaking, evolutions of NCA are 
not uniquely determined by initial configurations. A behaviour of NCA may be 
described by a state transition network (look at [4]). It is a graph, each of whose 
nodes represents some configuration. Directed arcs join the nodes to represent the 
transition between configurations. All nodes have out-degree one iff the cellular 
automaton is deterministic. Moreover, if the cellular automaton at hand is in fact 
nondeterministic and we consider the configurations in unlimited space U , then 
out-degrees of some nodes in the state transition network will be infinitely large. 

The difference between one and high dimensional CA has emerged from the 
solution of the predecessor existence problem (PEP) for CA. It is the problem of 
existence of a predecessor for the given configuration. S. Wolfram showed in article 
[8] that PEP is decidable for one dimensional CA, and T.Taku in article [10] showed 
that PEP is undecidable for some CA of dimension r where r — 2,3, ... . 

The computational equivalence of CA and Turing Machines is a well known fact 
(look at [2], [11]). The time complexity T(n) and the space complexity S(n) can be 
defined routinely for any CA A. 

This brings up the question: given an arbitrary predicate V , how does the min- 
imal complexity of r dimensional CA, computing V , depend on r? 

More precisely, if some predicate V is computable on CA of dimension r with 
time complexity T{n) > 0(n) is it possible to compute V substantially faster on 
CA of dimension r' > r ? 
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This is called TCD-problem. This problem is open for CA. 

A different situation arises with TCD-problem for NCA. For example, it is found 
that if some predicate V can be computed on NCA with time complexity T(n) — 
0(n a ) and space complexity 0(?W 2 ) then the increase of dimension by one unit 
allows to compute V in time 0(n a ^ 2 ). Moreover, given > 0, the time complexity 
0(n^) can be attained for the computation of such predicate if we increase the 
dimension of NCA to a suitable value. A similar result takes place also for faster 
increasing functions T(n). It means that multi-dimensional NCA will become the 
faster instrument for computation as the dimension increases. 

We proceed with the exact definitions. All constants are assumed to depend on 
the dimension r. 



2. The main definitions and results 



Let w — {c , . . . , Cfc} be an alphabet for possible states of cells. Let t takes the 
values from the set N = {0, 1, . . . }. 

An evolution in IT is a function of the form a : N x Z r — ► u>. A configuration 
is a function of the form : Z r — ► u. Any evolution a may be displayed as a 
sequence 

(1) a(°),a (1) ,-.. 

of configurations at the instants of time t = 0, 1, ... , where a^\i) = a(t, i). jth 
componenet of i e Z r will be denoted by ij. If i — (ii, i 2 , ■ ■ ■ , v) € Z r , we shall 
write a(t, i\, . . . , i r ) instead of a(t, i). 

The following notations are fixed for the cells comprising the neighborhood 
of the cell i: 

i(0) - i, 

i(l) = (n - l,i 2 , ■ ■ -,i r ), 
i(2) = (h + l,i 2 , ■ ■ . , i r ), 
i(3) = (ii,i 2 - 1,«3, • • - ,ir), 
i(4) = (i 1 ,i 2 + 1,23, • • -,i r ), 

i(2r) = (ii,i 2 , ...,i r + l). 
We put : a,j(t,i) — a(t,i(jj). 

We'll consider only such evolutions a that 3C : Vi : ||?'| > C a(t,i) = c , 
therefore all configurations will be finite objects. 

Let l r be a fixed computable one-to -one mapping l r : N — ► Z r such that 
||Z r (n)|| = 0(n 1//r ). This function represents an embedding of one dimensional 
space into r dimensional with the least norm ||Z r (n)||. 

Given an alphabet a, the set of all words over a is denoted by a*. If B = 
Cj 1 Cj 2 . . .Cj g G us* , then the initial configuration corresponding to the word B is 
defined by 

a (°)(i) = { ° 3k if ° " = k ~ S ' 

} c otherwise, 

we denote this configuration by a^P . 
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A nondeterministic cellular automaton of r dimensions is a function of the form 
A: ijJ x u> x • • • x lu — > 2". 

S v ' 

2r+l times 

Let b G o) 2r+1 , c G ^4(6). Any word of the form — > c is called a command of 
this automaton. This command is called trivial if A(b) = {c} and b has the form 
(c, tii, • • • , ^2r)- A set G of commands of A which contains all nontrivial commands 
is called a program of A. The behaviour of A is defined by its program. 

If V6 € w 2r+1 and A(b) consists of exactly one element, then in fact A has 
the form u 2r+1 — ► u, and we obtain the definition of a deterministic cellular 
automaton. 

We assume that A(co, . . . , c ) = c . This letter Co plays the role of blank, it is 
denoted by 0(w). 

An evolution of NCA A is a sequence a of the form (1) where Vt = 0, 1, . . . , \/i G 

(2) a(i + e .4(ao(M), ai(t,i), . . . a 2r (t, i)). 

It is obvious that in deterministic case depends on A, t, a^-°\ and in nonde- 
terministic case, in addition, on the choice of elements (2). 
The set of all values of a function / is denoted by Im/. 

Let the alphabet u> be divided into two nonintersecting parts: lu = w'Uw", where 
u' is the set of main letters, to" is the set of auxiliary letters, and let E C u>" be 
the set of end letters, where Ck-i,Ck € E. We denote u>'* by S, by succ(yl), E 
by Let for the evolution (1) of NCA A Im C w'U{c }. Let r(a) be the 

least value of t such that there exists one and only one letter c € E Him . This 
letter c is denoted by res(a, A) and is called the result of the operation of A on the 
initial configuration in evolution a. A configuration a r / a ) is called a resulting 
configuration for . 

In general terms, the result of the operation of A is defined uniquely only in 
the deterministic case. The set of all results in evolutions which begin with o^ ^ is 
denoted by yt[a( )]. 

A predicate P on the set S is an arbitrary subset of S. 

NCA A computes a predicate P iff \/B e S 

svlcc(A) e ^[a^], if B e P, 
succ(^) ^ ^[a^ } ], liB^P. 

It's obvious that such a predicate P is defined uniquely for A if it exists. We 
denote this predicate by P4. A cell i G Z r is called accessible in evolution a iff 
3i' < t(o) : a(t',i) ^ c . 

The diameter of the set of all accessible cells is denoted by D(a). 

Given B G P, the least value of r(a) from all evolutions a , where = 
a!§\ res(a,^4) = succ(yl) is denoted by tj^(B). 

Let Da(B) denotes the least value of D{a) from all evolutions a where = 
ves(a,A) = succ(^4), t{o) = ta{B). 
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The time complexity of NCA A is the function T A : N — > N defined by 
T A (n) = max{T A (B) \ B E P A , \B\ < n}, 
where \B\ denotes the length of B. 

The space complexity of A is the function S A : N — > N defined by 
S A (n) = max{D A (B) \ B e P A , \B\ < n}, 

Without loss of generality we may anticipate that A(ci, . . . ) = Cj for all Ci e E. 

Functions T A (n) and S A (n) can be very complicated, and it's convenient to use 
their best upper approximations T A , S A € K : T A {n) = 0(T A (n)), S A (n) = 
0(S A (n)) instead of them, where JC is the class of functions constructible in time 
(see below). Given the space Z r , a function / : N — ► N is called constructible 
in time on r-dimensional NCA, if there exists a constant c and NCA A with time 
complexity T(n) = c(f(n) + n 1 /*") such that for every word B e cj* , \B\ = n there 
exists the single resulting configuration for which has the form 

_ r C i,ifie{l,2,...,T(n)/2K, 
a T (t) = < 

Co, m the opposite case. 

It means that r-dimensional cube of side T(n)/2 can be isolated in time T(n), 
where 0{n x / r ) is the size of necessary domain for input word B. 

For example, the constructibility in time for n A and 2 2n is in fact proved in the 
section 4 (group Gl), for the functions n a ,q an , q, a e Q and for their combinations 
with additions, multiplications and superpositions the constructibility in time may 
be proved along similar lines. 

A pair of functions {T A ,S A ) is called a complexity of NCA A. Thus in what 
follows T A , S A (or, simply T, S with or without indices) will be constructible in 
time. 

The class of predicates P, computable on NCA of dimension r with complexity 
(T, S) is denoted by NC(r, T, S). 

We'll write T x < 0(T) instead of the following: VC > 3N Vn > N Ti(n) < 
CT(n). 

r-dimensional NCA, computing predicate P with complexity (T, S) is called a 
fastest NCA if P can not be computed on r-dimensional NCA in time T\ < 0(T). 

Let T^p denote the time complexity of a fastest r-dimensional NCA computing 
predicate P. 

Here are the main results of this paper. 

Theorem 1. 

NC(r,T,S) C NC(r + 1, VT + S, VT + S). 

Theorem 2. 

NC(r,T,S) CNC(r, Tx.T!), 

where Ti = 

Theorem 3. Let A be a fastest NCA of dimension r. Then 
(3) T A (n) = 0{S A (n)). 

Theorem 4. 

1) . T r _ liP = 0((T r , P ) 1 + 2 /'-), 

2) . T r+1 ,p = 0((T r , P ) 1 -'-/( f -+ 1 ) 2 ). 
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Now we shall give the outline of the following sections. All these results are 
based on two main methods of speeding up computations: The method of direct 
simulation in r + 1 dimensional space (section 3, Proof of Theorem 1) and the 
method of optimization of NCA in the same space (section 4, Proof of Theorem 
2). Theorem 3 will be simply derived from Theorem 2. Point 1) of Theorem 4 
will be proved by the method of simulation in r — 1 dimensional space (method of 
evolvents, section 5). Point 2) of Theorem 4 will be proved in two steps: reduction 
of space complexity in r + 1 space and the following optimization. 

Note that the both two method of speeding up require nondeterminism. 

3. Simulation in r + 1-dimensional space. Direct method 

Proof of Theorem 1 . 

Let P = P_a £ NC(r, T, S), A be a cellular automaton of dimension r with 
alphabet uo and complexity (T, 5), 0(w) = c . 

In this section we'll present the direct method of speeding up: we shall construct 
the nondeterministic cellular automaton NCA1 of dimension r + 1, which simulates 



The rough idea is that we expand the alphabet of the cellular automaton A at 
hand and use r + 1st dimension to code H state transitions of A into one big r + 1 
dimensional state transition of the new automaton NCAl simulating A, where H = 
0(VT + S). The single obstacle which will remain is that the initial configuration 
for NCAl is not but the ascending map, corresponding to input word B (the 
definition is in the next section). This obstacle will be overcome in the last part of 
this section. 

Definition A port is a list p of the form: 

(4) p = (Mark(p),con(p),envi(p),env 2 (p), . . . , env 2r+2 (p)), 

where A, B are the special new letters, Mark(p) £ {A, D} , the other members of list 
(4) are arbitrary letters from u>, and the functions con(p), env^ (p), j = 1, . . . , 2r + 2 
are defined by equality (4). 



Let <x>o be the set of all ports with the exception of (D, cq, . . . , cq). Then the 
alphabet of NCAl is ui\ = ujo U {®,b}, where 0,6 are new letters, and let 0(u>i) = 
(A, c , . . . ,co). 

Commands of NCAl 

Since the dimension of NCAl is r + 1, the commands of it have the following 
form: 

(P0,Pl, ■ ■ ■ ,P2r+2) ► M- 

Let allp q e uj\, q = 0, 1, . . . , 2r+2. The set {p} consists of all elements p £ lo\ to 
be described below. Let s+ = {0, 2, 3, 4, . . . , 2r + 2}, s~ = {0, 1, 3, 4, 5, . . . , 2r + 2}, 



A in time 0(T 



a / 2 + S). 



Alphabet of NCAl 




j + 1, for j odd, 
j — 1, for j even. 



We'll consider separately five different cases. 
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Case 1. Let the following conditions be satisfied: 



Vj e s + pj e wo, Mark(pj) = A, 
Vj £s + - {0} confjjj) = envj(p ), con(p ) = env k (pj), 



where k = k(j) and 

con(p 2 ) G -4(con(p ), con(p 3 ), . . . , con(p 2r+2 )). 

In this case p is such a port that Mark(p) = D, or p is po- 
Case 2. Let the following conditions be satisfied: 

Vj e s~ e wo, Mark(pj) =Dor pj = 0(a>i), 
Vj e s~ - {0} con(pj-) = eny,(p ), con(p ) = env fc (pj), 
where k = k(j) and 

con(pi) e -4.(con(p ), con(p 3 ), . . . , con(p 2r+2 )). 

In this case p is such a port that Mark(p) = A, or p is po- 

Case 3. Let p\ = 6, Mark(po) = D or p 2 = b. 

Then p is po or p is obtained from p by the following redefinition: we put 



Case 4. If 3j G s + n s : pj = b, we put p = b. 

Case 5. In all other cases we put p = 0. 

Note that NCA1 is nondeterministic even though .4 may be a CA of deterministic 
type, because in the cases 1 and 2 the choice of p is not uniquely defined. 

Let t be time step, H be some positive integer. We suppose that H is arbitrary 
till Lemma 3. 

Some peculiar configurations of NCA1 are called H, i-maps. Maps will be of two 
sorts. 

Definition. An Ascending map of high H at time step t (A, H, t -map ) is such 
a configuration a® for NCA1 that Vi G Z r+1 

1° If ii £ {-1,H}, then oW(i) G ujq . If a(t,i) 7^ 0(wi), then 1 < i± < H, and 
a(t, -1, i 2 , . . . , v+i) = a(t, *2, • ■ • ir+i) = b; 

2° Hi: < h < H- 1, then 

con(a(i, ii + 1, i 2 , i 3 , . . . , v+i)) = -4(con(a (t, i)), con(a 3 (i, i)), . . . , con(a 2r+2 (i, i))), 



then Vj : 1 < j < r+1 3e e {I, —1} con(a(t, i\, . . . , ij+e, . . . , V+i)) = cnvfc(a(i, i)), 
if the two parts of this equality exist. 




3° If 
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Definition. A Descending map of high H at time step t (D, H, i-map) is such 
a configuration for NCA1, that Vi G Z r+1 
1° See above. 
2° If < ii < H, then 

con(a(i,ii - 1,4 2 ,»3 ) • • •, V+i)) = -4(con(ao(t,i)),con(a 3 (t,i)), . . . , con(a 2r+2 (t, i)). 

3° See above. 

The following proposition relates evolutions of A to those of NCAl. 
Lemma 1. of ' a< T ',... «s 

evolution of A iff there exists an evolution 

(5) a<°>,...a« ... 

o/ NCAl, where for every r = 0, 1, . . . i/ie following condition C T is fulfilled: 

C T . 3t = t(r) eN3q = q(r) : 0<q<H-iy~ieZ r con(a(i, q, h, i 2 , ■ ■ ■ , v)) = 
a(r, i), w/iere /or i even; r = — 1) + q, and is A,H,t-map, for t odd: 
t = (t + 1)(H - 1) - q and a (t) is D, H, t-map. 

Proof. 

1. Necessity. Let a be an evolution of A. An evolution a of NCAl is called 
r-evolution if the condition C T is fulfilled. At first let us prove by induction on r 
that for any t there exists r-evolution. Basis: r = follows from the definition of 
A, H, i-map. Step: follows from the definition of NCAl. □ 

Now we introduce the order -< on N x {0, 1, . . . , H — 1} by the following: (ti , q\) -< 
(t 2 ,q2) iff t\ < t 2 or (t\ = t 2 - even and qi -< q 2 ) or (ti = t 2 - odd and q 2 ~< qi). 
Note that if n < r 2 , then (t(ri) , q(ri)) ^ (t(r 2 ) , q(r 2 )) . Consequently, in view of 
the definition of r-evolution, if a is Ti-cvolution, d is T2-evolution, then for every 
pair (t, q) ^ (t(ri), q(ri)) we have Vz e Z r a(t, q, i) — d(t, q, i). Thus we obtain that 
there exists the evolution a of NCAl which coincides with every r-evolution on all 
lists (t,q,i), where (t,q) ri (t( T ) , l( T )) ■ Necessity is proved. 

2. Sufficiency. Follows from the definition of NCAl. Lemma 1 is proved. 

Lemma 2. If , . . . a^> , . . . is evolution of NCAl and is H,0 -map, then 

aW is H, t-map iff 

(6) VieZ r+1 : 0<h<H-l o(t,i)^0. 
Proof. 

Induction on t. Basis: t = 0. Follows from the condition. Step follows from the 
definition of NCAl and those of H, t-map. Lemma 2 is proved. 

Let [x] denote integral part of x e M. 

Lemma 3. Let some predicate P be computed on NCA A with complexity (T,S), 
\H — [yjT{n)\\ < yjT(n)/2 , B e E, a' ^ be the initial configuration of A, corre- 
sponding to B,t = 2[ V /T(n)], n=\B\. 

If B £ P then for some evolution of NCAl of the form (5) the condition (6) is 
fulfilled and 

(7) 3i e Z r+1 : con(o(t, i)) = succ(yl) 

and if A ^ P then for any evolution of NCAl of the form (5) with the condition 
(6) the following property takes place 

(8) Vi : con(a(t, i)) ^ succ(*4). 
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Proof. 

Follows immediately from the definition of computation on NCA, Lemma 1 and 
Lemma 2. Lemma 3 is proved. 

4. Auxiliary automata 

To finish the proof of Theorem 1 we must construct two auxiliary NCA: 

1) . NCA2 - which begins to operate from d^ - initial configuration, corre- 
sponding to B e £ in Z r+1 , so that for some evolution a from Lemma 1 and for 
some H = a(°), \H - [y/T{nj}\ < y/Tfnj/2, n=\B\. 

2) CA3 - deterministic CA, which begins to operate from arbitrary chosen H, t- 
map from (5) and has the result only if (6) is fulfilled, and in this case 

f succ(CA3) e CA3[a,W], if (7), 
\ succ(CA3) i CA3[a { % if (8). 

If now we combine the sets of commands of NCA1 , NCA2 and CA3, and put 
succ(£>) = succ(CM3), E(B) = E(CA3), then in view of Lemma 3 the resulting 
NCA B will satisfy the conditions of Theorem 1 . 

Now let us describe NCA2 and NCA3. Real work of NCA2 is such that at 
first the group Gl will operate, and after that, groups G2 and G3 will operate 
simultaneously. 

NCA2. 

The program contains 3 groups of commands: 

Gl : Realization of function H = 0(^T A (n) + SU(n)). 

G2: Isolation of area for modeling. 

G3: Construction of A — H — 0-map. 

If , for example, T A (n) = 2 2 ™, S A (n) 
the following forms: 



(a, . . . 


,2/,?) — ► Ci, 


(c ,c , . . 


■Ci) 


— >d, 




!)•■•) > Cj , 


(c i; . . . 


,co) 




(Cj, . . . , 


Cj,?) — ► c*, 




...) 






Cj, ?) > Cj, 


{pi, ■ ■ ■ , 7^ C s 


.4) 




(d, c , . . 


. . ,Cj) — > d', 


(y,.-- 




— ► e, 



where s take all values from {1,2,...}, ? means an arbitrary letter, y € 
{co, rf, d'} , Cj, c*, c^ are new different special letters for any Cj G lo. 

Group G2 consists of all commands of the following forms: 

(c Q ,d',...) — *d, (e, ...) — > c , 

(d, d',...) — yd' (zi, c , z 2 , • • • co) — > c , 

(d',^d', ...) — > d, (zi,c ,c , ■ • • ,co) — yb, 

(...,&,...) — > b, (c ,c ,e, ...) — > 6, 



where 21,22 € {d, d'}. 
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Group G3 consists of all commands of the following forms: 



(cj,b,. 
{co,hi, . 
(hi, . 
(hi,, 
(h 1 ., • • 



(r 



3 ' ' 



• h» 



hj, 
h s , 
hi, 



where i, j, s take all values from {1, 2, . . . }, hi are special new letters, corresponding 
to each i, and Ti takes values from the set of all ports p such that con(p) — Ci. 

For Tj((n) = n 4 , SU(n) < 0(VT) the group Gl consists of all commands of the 
following forms: 



(co, ■ 



(ci, ...,c ,^ c°) 
(c 4 + ,---,^c ) 
(Ci, . . . ,c i , ?) 
(C*, ■■■,¥= Co) 
(c*, . .. ,c ) 
(Ci,-- - t^c+,cT) 
(ci, . . . Cj , 7^ c s ) 



/ 



(cr....^^,?) 




(c-,..., c ;,?) 




(c'i,...,^cj) 




(c'i,---,cj) 




(c+ ...,co) 




(c?,...) 




(Ci,...C?) 




(*,■■■<£) 


— ► e, 


(rf, . . . Ci) 





where z e {e?, e?'}. NCA2 is described. 
CA3. 

PutE{CA3) = {cf k ,d k _ 1 }. 



(r Cp ,^0,6,^) 

no 



i = 0,l,...,2r, P e{M-i}, 



2r+l 



(...,0,.. 



where every letter of NCA1, NCA2 may occur in " ..." . 
Theorem 1 is proved. 

The general form of "successful" operation of resulting cellular automaton is 
shown on Figure 1. 

The technique evolved allows to derive the following amplification of Theorem 

1. 
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Corollary 1. If r 2 > r\, then 

NC(n,T, S) C NC(r 2 ,VT+S, (S^T 1 ' 2 ) 1 '^ + S). 

Proof 

Let A be a nondeterministic cellular automaton of dimension ri computing 
predicate P4 with complexity (T,S), S(n) < 0(T(n)). Nondeterministic cellu- 
lar automaton of dimension r 2 simulating A with complexity (T^Si), T\ = 
0(VT + S), Si = 0((S ,ri T 1 / 2 ) 1/r2 + S) can be constructed by a modification of 
method from the proof of Theorem 1 . 

The area for simulation can be organized so that its diameter will be 

0((5 ri T 1 / 2 ) 1 / r2 + S) 

For example, for n = 1, r 2 = 2 the area from the Figure 1 can be constructed 
as a spiral so that the operation of B will have the form, presented on the Figure 2. 

Here the work of B in straight strips is the same as in the proof of Theorem 1, 
and in corner squares B only verifies the coincidence of the words written on two 
sides at each time step, like on Figure 3. It is not difficult to understand how B 
must be arranged, and I omit awkward details. □ 

Corollary 2. // InT/ \nS = 0(1), then for all r = 1,2, .. . 

NC (r, T, S) C NC ( [In T/ In S] + [log 2 (In T/ In S)] + 2 + r, S, S) . 

Proof 

Lemma 4. Let A be NCA, computing P4 with complexity (T, S*), S < 0(T). For 
every k = 0, 1, . . . there exists a sequence 

Ao, A\, A 2 , ■ ■ ■ ,Ak, 

where A = Ao, for all i = 0, 1, . . . Ai is NCA of dimension ri and complexity 
(Ti, Si), computing Pj,, where for i > 

n < r,_i + pnTi/21nS] + l, 
(9) = 0(^^77+5), 

Si - O(S). 

Proof 

Induction on k. Basis is evident. Step. Let k > 0. Applying the inductive 
hypothesis we obtain the sequence Ao, A\, A 2 , ■ ■ ■ ,Ak-i with the conditions (9). 
When T fe _i = 0(S k -i), we put A k = A k -i- When T fe _i > 0(S k -i), Corollary 1 
with Tfc_i playing the role of T yields the required NCA A k . Lemma 4 is proved. 

Let h be the least number such that Th = O(S). In view of Ti = 0(-v/Tj_i + S) 
we conclude that h < [log 2 (lnT/lnS')] + 1. Then it follows from (9) that 

r h < r + h + 1 + L/2 + L/4 + • • • + L/2 h <r + h+l + L, 

where L — [lnT/lnS*]. Corollary 2 is proved. 
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4. The optimization of NCA. Fastest NCA 

Proof of Theorem 2 

Proposition 1. Given a function a(n) : N — > N, constructible in time , every r- 
dimensional NCA A with complexity (T, S) can be simulated in IT with complexity 
{C 1 T/a + C 2 Sa 1 / r , CaSa 1 ^), C U C 2 ,C 3 depend on r. 

Proof. 

As above, we shall consider the case r = 1 in more detail. NCA A', simulating 
A has the form 

A' =Ai *A 2 , 

where A\ constructs the area for simulating and A 2 simulates A. 
Given an input word B , the area for simulating has a form 

(10) ...c o bB0 k {b0 p ) a b lCa ... , 

where p = k + \S\. The configuration (10) can be simply constructed in a time 
0(aS) by A\. The configuration (10) is an input one for A 2 . We suppose that 
T/a = q{n) e N and let an evolution a^^a^, . . . of A be divided into a 

sequential segments: 

A i: a« ...,a(« +1 \ 
A 2 : a^ +1 \...,a^ +1 \ 

A a : a^ 1 ' ad, 

such that input configuration of Aj+i'- in(j + 1) and output configuration of A r : 
out(j) are identical, j — 1, . . . ,a — 1. 

A 2 simultaneously simulates all possible segments Ai,A 2 ,...,A Q at the sites 
of sequential occurrences of the words: B0 k ,b0 p , b0 p , ... in (10). For every Aj A 2 
stores in(j) and checks the equality 

(11) in(j + 1) = out(i) 

for all j = 1, . . . , a — 1. A 2 places a special mark of "flow" in a cell where the 
violating of equality (11) is detected. At last the letter b\ moves to left through all 
domains where the equality (11) has been verified and if b\ has not met a mark of 
"flow" before Co, then A 2 achieves "success" . All these actions can be fulfilled simul- 
taneously because A 2 is nondeterministic and we obtain that A achieves "success" 
in time 0(T) iff A 2 achieves "success" in time 0(T/a) beginning with (10). 



COMPUTATIONS ON NONDETERMINISTIC CELLULAR AUTOMATA 13 



A program of A2 has the following form: 

(c,?,?)— >(c,c) ) 
(0,?,?)— ►(a,*), 
((ci,c 2 ), (ei,e 2 ), (91,92)) — > (ci,^(c 2 ,e 2 ,32)), 
((*,y),?,?) — 
((z+y+),?,6) —►(*+()), 
(6, (s + ,tf + ), ?)—>»+, 
((* + ,0),y+?)— >(*+ 
(!,+,?,?)—> 6, 
(( 2 ; + ,0),(y + ,z + -,?)^( a; + ,z + -), 
((x+,y+-),?,(y+,0))^( a ; + ,0), 

(04 ,2/r),04,y 2 ~ )>04 )) — ► fr> 

where h = /, if x\ ^ y 1} or x 2 ^ y 2 , or x 3 ^ y 3 , 
else h = s', 
((x + ,y + ),u,v) — ► /, where e 
((a; + , ? y-), S ',fe)^, 
(6, s",c ) — ► s, 

((x + , y~) or b or s' or s", ?, s) ► s, 

((c,c),x,?) — ► /, 
((c,c),?,ar) — ► /, 
(s,co,?) — ► !, 

where succ(*4 2 ) = !, / denotes "flow". The case r = 1 is considered. In the cases 
r > 1 the input configuration Do for „4 2 consists of a copies of cube {1, . . . , S(n)} r , 
disposed sequentially along a spiral, so that the size of D does not exceed 0(Sa 1 / r ). 
Evident changes must be done in the definition of A 2 . Proposition 1 is proved. Let 
us turn to the proof of Theorem 2. 

Proposition 1 makes it possible to achieve the time complexity 

T(a) = d- + C 2 Sa 1/r . 
a 

With the aim of finding the minimum of T(a) we consider the equation: 

T'(a) =-d^ + CzSa-r^/r = 
or 

which yields 

(12) a min = C(T/S)^tt 

for some constant C = C(r). The function (12) is constructiblc in time as T and 
S. Taking this value of a for A' from Proposition 1, we obtain NCA A' simulating 

1 r 

A with time and space complexity 0(T~ S~ ). This proves Theorem 2. 
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Proof of Theorem 3. 

The assumption that T > 0(S) for some fastest NCA contradicts to Theorem 
2. Theorem 3 is proved. 

Hence, while on the subject of fastest NCA we may talk only about their (time) 
complexity T, because in view of Theorem 3S = 0(T). 

5. Complex method of simulation. Method of evolvents 

Proof of Theorem 4- 

Proposition 2. If r > 1, then NC(r, T(n), S(n)) C NC(r- 1, T(n)S(n), S^(n)). 
Proof. 

Given r-dimensional NCA A r with complexity (T, S) , computing predicate P, 
we shall define NCA ^4 r _i of dimension r — 1 with complexity (T r _i, SV-i), such 
that each time step in the evolution of A will be simulated by A r -\ in 0(5(n)) 
time steps. 

The support of configuration a^> is the set = Z r |a(t, i) 7^ c }. Without 
loss of generality we can assume that all supports Sr for A r are r-dimensional cube 
B r , B = {1, . . . , iJ}, where = S(|B|) depends on the input word B. All supports 
s£-i of corresponding evolution of A r -i are contained in r — 1-dimensional cube 
Bl' 1 , B 1 = where if p = min{pe N | pH > ff^r}then#i = (p+T)H. 

If Vf, denotes a section of B r by hyperplane i r = x, x = 1, . . . , H, then there 
exist Vi, . . . , Vh C Bl^ 1 and fixed isomorphisms V x — > V4 such that 

H 

Definition. The set (J V4 is called an evolvent of £T. An evolution of A r -i 

in V4 will simulate an evolution of A r in V' x . 

We shall describe A r -\ in detail only for r = 2 because for r > 2 A r -\ can be 
constructed along similar lines. To facilitate further notations we need to introduce 
some auxiliary notions. Let A and B be one dimensional CA with alphabets uj\ 
and LO2, determined by programs and Ub respectively. 

Definition. 

1) . A composition of A and B is a cellular automaton, denoted by A * B which 
is determined by program U fig. 

2) . A direct product of A and B is a cellular automaton denoted by A x B with 
alphabet oj° = (w\ x C02) U {0(o>i)} and program II, which consists of all commands 
of the form 

((ui,u 2 ), {vi,v 2 ), (wi,w 2 )) — > (zi,z 2 ), 

where z x = A(u 1 ,v 1 ,w 1 ), z 2 = B(u 2 , v 2 ,w 2 ) , 0(lu°) = 0(o>i). 

3) . In addition, let G be some set of words of the form: 

((ai,a 2 ), (0i, fa), (71,72)) — > (<5i,^), 



where ctj, 0j,jj,5j e ujj, j = 1, 2. 
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NCA V , determined by the program II U G is denoted by A xg B. It may be 
called a semidirect product. 

4) . Let IT^ be the set of all commands of the form (x,y, z) — ► u, where 
(x, z, y) — ► u contains in II4. Then A* denotes NCA with program 11^. 

5) . At last let Af u denote a standard automaton in alphabet u> : Af u (x, y, z) = x. 

Definition of A r -\. 
Put: 

Ar-l = B r -i * M, 

where M. is a marker of initial evolvent as and B r -\ simulates A r -\ in evolvent. 

Let uu = {co, c\, . . . , Ck} be alphabet of A r , uo + = {c^, c+, . . . , c^} and oj' = 
{cq, c[, . . . , c' k } be new different copies of u; 0, b be new letters. Put a ~ uj U uj + U 
Lu'U{0,b}. 

An auxiliary CA W in alphabet a is defined by the following program II: 

(d,?,x) — +0, (cj,?,0)— »0, 
(0,*,?)— (0,y,?)— y, 
(& )Ci ,?)— c+ ( C + ?,?)—>&, 

where x takes all values from {0, 6}, y - from lo + U a/, i = 0, . . . , k. 

Let uji = (ax ax er) U {c }, G be the list of all commands of the following form: 

{(U 1 ,U2,U S ),(V 1 ,V2,V 3 ),(W 1 ,W 2 ,W S )) ► (Z,Z,Z), 

where z = A r (u2 1 v 2 ,W2,u 3l u 1 ). 
Put 

B r _! = (Wx^)x G r, 

SUCc(S r _i) = (Cfe,Cfc,Cfc), 

£(B r _ 1 ) = {(c,c,c) | cetfpl,.)}, 
0(wi) = c , 

where Ck = succ„4 r , Co = 0(u>). 

Given an input word B S to* for *4 r , let d(r, s) be the initial configuration of A r , 
corresponding to B, we define corresponding input word as — h hi . . . h H 2, H = 
S Ar (\B\) for B r -i by the following: 



hi = < 



(b,b,b), if i = ( mod if), 
(c , c , c ), if i = r ( mod if), q < r < H, 
d(r, s), if i = si? + r, 1 < r < q. 



The following Lemma 5 may be deduced immediately from the definition (13). 
Lemma 5. VB e P Ar 7» p _i(ob) < 2Hr Ar (B). 

Now to finish the proof of Proposition 2 it is sufficient to construct CA M. 
transforming B to in time 0(H 2 ) = 0(H T ^ T ). It may be done as in the 
proof of Theorem 1 (look at NCA 2: Gl and G2). Note that here we have V x = 
{{x-l)H + l,...,xH}. 
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The case r — 2 is considered. 

Let now r > 2. This case differs in that the domains V x for sections i r = 
x, x — 0, 1, . . . , H, H = O(S) are disposed sequentially along a spiral so that 
S r _i = O(S^r), T r _i = 0(TS). Such initial evolvent a R C T" 1 can be isolated 
from the space and marked out according to the direction of laying of all V x into 
evolvent in time 0{S). 

Proposition 2 is proved. 

Note that analogous proposition takes place also for deterministic CA. 
Now we turn to the proof of Theorem 4. 

1) . Given a fastest r-dimensional NCA A with complexity T — T^p, P = 
P4, Proposition 2 yields r — 1-dimcnsional NCA A' simulating A with complexity 
(T 2 ,T^t). Then by Theorem 2 we obtain r — 1-dimcnsional NCA simulating A' 
with complexity 0{T 2 ' r {T^)^) = 0{T 1+2 / r ). Point 1) is proved. 

2) . We can suppose that A acts in r-dimensional cube B r , B = {1,...,S}. 
r + 1-dimcnsional NCA simulating A in time 0(T 1 ~ r /( r+1 ) ) will be constructed in 
two steps. 

Step 1. Simulation of A in Z r+1 with complexity (T, T^tt). 
Step 2. Applying of Theorem 2 to NCA obtained in Step 1. 
Step 1. 

Definition. A set A consisting of inclusions of the form 

Lj : Y j — > W, j e N, 

is called a constructible set of inclusions if for some q < r all Y 3 ■ c 7L q , and for 
some constants c, ci , c 2 the following three conditions are satisfied. 

a) . Vj g N, x,yeYj 

(14) p(Lj(x),Lj(y)) < cp{x,y), 

where p denotes standard metric in Z 9 or 7L r . 

b) . Diam(Imi :) ) < c\\Lj\^l r , where \L\ denotes the number of elements in L. 

c) . All domains ImLj can be marked out by NCA in time C2\Lj\ 1 / r so that every 
cell z = Lj(x) G ImLj will be marked by a label m(z) which points to the dispo- 
sition of all such cells Lj(y) that p(x,y) = 1 with respect to z. (Every label m(z) 
has the form (g ,gi, ■ ■ ■ ,g~2q) where g p = L (x{p)) - z G Z r , p = 0, 1, . . . , 2q; x(p) 
is defined in the section 1. In view of inequality (14) the required number of all 
such labels m(z) does not depend on j.) 

To fulfill Step 1 it is sufficient to prove the following 

Proposition 3. For every r there exists a constructible set of inclusions L r s : 
B r — > Z r+i^ B = {1, . . . , S}, S = 1,2, 

Really, with such inclusions L r s in view of the point b) we can simulate A with 
complexity (T,0(S*r)) in Z r+1 . 

Lemma 6. Given r < 1 there exists a constructible set of inclusions 

R r s : Y s ^Z r +\ 5=1,2,..., 



where Y s = B{ x B G Z r+1 , B x = {!,..., Si}, Si = [S^]. 
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Proof. 

A parallelepiped Yg in Z r+1 can be conceived of as a thread of length S and 
S~^~ thick. This thread can be rolled up into a ball which has the shape of cube of 
side 4|F5|~ in Z r+1 , because < \Y S \~ = S* 7 ^ . The construction of NCA 

performing the required marking of ImYs is evident. □ 

Proof of Proposition 3. 

Induction on r. Basis r = 1. Proposition 3 follows from Lemma 6. Step: r > 1. 
It follows from the inductive hypothesis that there exists a constructible set of 
inclusions: 

Vf 1 : B 7 - 1 — ►Bf, S = 1,2,.... 

Let R r s be inclusions from Lemma 6. Then we obtain a constructible set of 
inclusions U' s : B r — > Z r+1 , defined by the following: 

L r s (x,y) = K r s (L r s -\x),y), 

where x <E Z r_1 , t/eZ. 
□ 

Step 1 is fulfilled. 
Step 2. 

Applying Theorem 3 we obtain a simulation of A in Z r+1 with complexity 
Theorem 4 is proved. 

Thus, if T4 > S 1 ^, then we have a variety of ways to accelerate the computations 
of Pa in r + 1-dimensional space, for example: 

1) . By Theorem 1 (or by Corollary 1). 

2) . Sequential applications of Theorem 2 and point 2) of Theorem 4. 

The second way gives the better acceleration for T4 3> SU- But if T is small, 
for example T = S 2 , then Theorem 1 gives the more strong estimate for T r+i ,p 
because 1 < ^±2(1-^-1^). 

Remark. For any predicate P computible in time T = 0(n a ), a = const on 
d- dimensional NCA and for any [3 > there exists such a number r > d that 
T r ,p = 0(r/). 

Proof. 

Applying Theorem 4, point 2), we obtain that for r > d 

r 

T r , P (n)<c(R)exp[MT d . P (n)) J] (1 - m/(m + l) 2 )]. 

m— d+1 

We have: 

r r 

In J] (1 - m/((m+ l) 2 ) ^2 m/(m+l) 2 — > -00 (r — ► +00). 

m— d+1 m— d+1 

Thus, f[ (l-m/(m+l) 2 ) — > (r — > 00). □ 

m— d+1 
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6. Discussion 

We see that if cellular automata are nondeterministic, then increase of dimension 
leads to the substantial acceleration of computations. In addition, programming on 
nondeterministic CA is more simple, then on ordinary CA. Thus, if it is possible, the 
realization of many dimensional nondeterministic cellular automata by a physical 
device would be of great practical consequence. From the other side, TCD- problem 
for deterministic CA remains unsolved. Let C(r,T,S) be the class of predicates, 
computable on CA in time T and space S. TCD-problem for CA is as follows: 

Given r,S,T, is there an increasing function f(n) such that 

C(r,T,S)CC(r+ l,^+5,5i) 

for some S± ? 

At last note that the hypothesis T = 0(f(S)) for the fastest deterministic CA 
is open for question for every function f(n) > n growing slowly in comparison with 
exponential. 
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